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ABSTRACT 

We analýze the correlation function of a meson with one heavy and one light quark in 
inverse powers of the heavy quark mass mg using a succession of Foldy-Wouthuysen-type 
transformations prior to radiative corrections. We evaluate the correlator to order itlq 
^ . in a random and dilute gas of instantons, using the planar approximation. We show, in 

leading order in the density, that the heavy quark mass is shifted to the order vplq and 
that the induced interaction between the heavy and light quarks is attractive. We also 
find it to be an order of magnitude smaller than the 't Hooft interaction between two 
light quarks. The shift in the heavy quark mass is related to the perimeter law of large 
Wilson loops. The relevance of these results for generál hadronic correlators with heavy 
quarks is discussed. 
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The instanton model for the QCD ground statě [|lj offers an interesting framework for 
the discussion of soft physics from first principles. While instantons do not confine color 
in pure Yang-Mills theory, they provide a simple mechanism for the spontaneous breaking 
of chiral symmetry in QCD. Extensive calculations using instantons both analytically and 
numerically suggest that the instantons in the vacuum form a dilute gas 0. This may 
be understood by noting that the instanton-antiinstanton interaction is not dipóle like at 
large distances but screened by the fermions over distances on the order of 1/2 fm. A 
dilute systém of instantons behaves as a free gas. As a result, the singlet rj' is heavy and 
the topological susceptibility is small. 

Hadrons with a heavy quark exhibit a new type of symmetry: invariance under spin flip 
of the heavy quark. This symmetry has been used in combination with chiral symmetry 
(Heavy Quark Effective Theory) [[| to understand the structure and decay properties of 
heavy-light systems. The main observation is that if the mass of one quark is taken to 
infmity and factorized, the physics of the remaining systém is soft and constrained by 
chiral dynamics. The analysis of the soft physics in heavy-light systems can be assessed 
quantitatively using QCD lattice simulations, or qualitatively using QCD inspired models. 

In this letter, we investigate the effects of a dilute gas of instantons on mesons com- 
posed of one heavy and one light quark. We show that the heavy meson correlator may 
be systematically investigated in inverse powers of the heavy quark mass using standard 
Foldy-Wouthuysen transformations @] and before the inclusion of radiative corrections. 
We analýze the heavy meson correlator in the planar approximation and for a dilute gas 
of instantons. To leading order in the density, the instantons shift the heavy quark mass 
as first noted by Diakonov, Petrov and Pobylitsa ||, and induce an attractive scalar-color 
interaction between the heavy and light quarks that is an order of magnitude weaker than 
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the induced 't Hooft interaction between two light quarks . The shift in the heavy quark 
mass is related to the perimeter law of the Wilson loop and reminiscent of the classical 
Coulomb energy. The instanton gas screens but does not confine. Our construction is im- 
mediately amenable to heavy baryons. The instanton size p is fixed to 1/3 fm throughout 
and the instanton density is fixed to 1 fm -4 . 

1. Consider the following meson correlator in Minkowski space 

C T {x, x 1 ) = (0|T {qY^{x)^Vq{x')) |0) (1) 

with T = (1,7) x (1,T) an arbitrary spin-flavor channel. In the regime where the mass 
rriQ of the field ip is heavy, the correlator (|l[) may be analysed systematically in 1/niQ 
by redefining the field ip using the Foldy-Wouthuysen (FW) construction 0. In the bare 
version and to leading order in l/mg, the transformation reads 

i/)(x) ~ e - Í7om< 3 ť e - ířT,M[v °' V!1/4m Q e "^ /2m « Q{x) (2) 

where V = d — iA and with g (the gauge coupling) set to one. The first transformation 
rescales the momenta, the second eliminates the odd parts, and the third removes the mass 
term. The successive transformations in (2) are vector-like, unitary and gauge-covariant. 
Thus, they are norm-preserving and anomaly-free. In terms of (0) the QCD part of the 
action for the heavy field ip becomes 

C^Qí^Q + q[-^- q -^Q (3) 

with B % = — ze lJ ' fe [V\ V J ']. Equation @ has the expected FW form to order mg 1 . Under 
(§) the heavy meson source shifts 

qTi) ~ qYe-^ m ^ [ 1 - ] Q (4) 
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As a result, the correlator ([!]) takés the generic form 



C r ~ (Olgre-^^-^QQrglO) + {0|gr e -^° m « (í+í ' ) 



QQ !Ť V 



2m c 



r g |o> (5) 



This construction can be extended to arbitrary orders in l/mg. Since the FW transforma- 
tion preserves manifest gauge invariance order by order in the heavy quark mass, the issue 
of renormalizability holds order by order in 1/tuq. Our construction refers to the bare 
quark mass rriQ. We are assuming that the expansion is not upset by renormalization. 

We note that the second term in ([5]) induces mixing between the particle and the 
antiparticle content of the correlator. This mixing follows from the redefinition of the 
sources and not the action. As we are interested in the case where rriQ is large ( non- 
relativistic limit ) it is convenient to take the projected sources T± = r(l ±7°)/2 instead 
of the unprojected sources T. Since r ± 7 • Vr± = the mixing part in (|) drops. Thus, 
in Euclidean space 

Ct{x,x') ~ -e^ T - T '\0\Tr(r±S Q (x,x')r±S(x' \x)) |0> (6) 
where S is the propagátor of the light quark, and Sq is the heavy quark propagátor, 

with = 74/ÍV4 being the free part. The order itiq has been considered by many 
authors Terms of order m^ 2 and higher can be sought along the same lineš. 

The present construction has the advantage of transforming the complete correlator 
and action into an effective correlator and an effective action with manifest power counting 
in 1/m.Q. The transformation does not rely on the commonly used momentům decompo- 
sition pq = txiqv + k, and as such does not suffer from the ambiguities associated with 
the redefinition of the heavy quark velocity v S. We have explicitly used the rest frame 



of the heavy quark, thus v = (1,0) in our case. The derivation can be easily generalized 
to arbitrary frames. 

2. We now proceed to evaluate (0) to order mg in a random gas of instantons. In 
generál, the correlator C ~ (S Cg) ^oo) receives contributions from both planar and non- 
planar graphs and thus is difficult to analýze. However, to leading order in 1/N C the 
planar graphs dominate |5], KJ. To this order, the r]' is massive (even in the chiral 
limit) and the instantons are screened. A dilute systém of instantons and antiinstantons 
with zero net topological charge behaves as a free gas. The systém screens but does not 
confine (see below). This notwithstanding, the planar graphs may be resummed using 
a Bethe-Salpeter-type equation. Following Pobylitsa fllOl , the inverse correlator @ after 
resummation reads 

c- 1 ~ s- 1 ®sl- 1 -J2((s-4 TY 1 ® (^oo - 4 sj) 7 "" 1 ) (8) 

The upper script T is short for transpose. In the planar approximation, the light propa- 
gátor satisfies the integrál equation 

s- 1 = s 1 -Y,((4l 1 -sy 1 ) (9) 

For massless quarks S^ 1 = +i$ , while the infinitely heavy quark propagátor satisfies 

s- 1 = s; 1 -J2((4ú-s oo y 1 ) (io) 

where 4 4 = 7 4 ^4 4 and S* = řy 4 <9 4 . The sum is over all instantons and antiinstantons and 
the averaging is over their position zj and color orientation Uj. Generically 

(vj éz ) /^+( J ^ 7 )~™/^ ( J + 7 ) en) 

where n = N/2V,iN c = n*/2N c . The last substitution follows from color averaging over 
two ř7/s for illustration. The diluteness factor of the gas is given by the dimensionless 
combination np 4 ~ 10~ 3 , with n ~ since n* ~ N c ||. 



In the dilute gas approximation, the nested integrál equations (||-|1Č]) may be iterated 
in n. To leading order 



<T 1 ~ s- L ® s. 



-1 CT,-1 



n ďzTrcUL]! ® [H]i + 1^1 



(12) 



with 



L ]/ <g> [i/]/ 



r,-l 
D 



|$o)($o| 



rr-1 
°0 



s- 1 



z 7 4 V 



5* s 



-1 
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Here |$o) is the normalized fermionic zero mode in the one-instanton background, and 
So = ($o|S|$o) ~ (240 MeV)™ 1 . We note that the effects of the planar resummation is to 
generate an effective mass for the zero mode in (|13|). This result is directly related to the 
induced constituent quark mass. Indeed, resumming the planar graph for the single quark 
propagátor shows that S^ 1 ~ Sq 1 + íy / n£, where £ is a space-spin-dependent kernel that 
satisfies the following gap-like equation JTO 



S = f*z Tr c (S ^* 0> <*° 1 & +1-7) (14) 

with S following from flíll) by taking the expectation value in the zero mode statě. The 
shift in the light quark mass is 



AM q ~ \pň < X-oolSlx+oo >~ y/n 



2J2 S 



420 MeV 



(15) 



Similarly for the heavy quark propagátor, we have ~ S# + in O, where is a 



t-dependent kernel satisfying 



9 = I ďz Tr c I S. 



-i 



i 7 4 V 



s* isz 1 + / - / 



4,i 



(16) 



In coordinate space, the heavy quark propagátor in the one instanton background reads 



< x\^—\0 >= 5(x) 9(t) Prf dsMxs-zj) 
iVai 



(17) 



with x s = (s,x). Inserting fll7|) into (0) and using the one-instanton configuration 

A;(x) = +% v x u ( 1 - -^iJ) (18) 



ar x z + 



yields for large times 



< X^elx+oo >~ (-327Tp 3 7 4 ) ( XC ° S ( 2 jíT^ ^) = (" 327r A 4 ^) (19) 

The corresponding shift in the heavy quark mass follows from the effective vertex (16) 

AM Q = 32vmp% ~ 16vmp 3 ~ 70 MeV (20) 

as first suggested by Diakonov, Petrov and Pobylitsa || . Note that the shift in the heavy 
quark mass is of order n as opposed to the shift in the light quark mass which is of order 
y/n, almost an order of magnitude smaller. This is simply the classical Coulomb energy 
of a heavy quark. A similar result has been recently obtained by Bigi, Shifman, Uraltsev 



and Vainshtein and also S. Narison |]TT|. Remarkably, their estimate using the infrared 
renormalon suggests AMq ~ 50 MeV. and AMq ~ 70 MeV respectively. 

3. The shift in the heavy quark mass to order mPq is related to the behavior of large 
Wilson loops in the random instanton gas. Indeed, the Wilson loop can be rewritten as 
follows 

(Tr c Pe^xL^- A ) = _^ T r c (W^i dx ' A{Tl) Pe^ dx ' A{Ll) Peh dx ' A{T2) Ve ^2 dxA{L2) ^ \ 

(21) 

where T 12 and L 12 are the paths shown in Fig. 1. A(T, L) = J2i Ai(T, L) is the sum of 
the instanton and antiinstanton gauge configurations projected onto the lineš T, L. For 
large Wilson loops T,L » 1/p, the heavy quarks (line integrals) no longer interact 
with each other. As a result, (^T|) factorizes into the trace of the product of four "heavy 
quark" propagators (Wilson lineš), each averaged over the ensemble of instantons and 



antiinstantons, and evaluated for large separations T and L respectively. The above 
discussion shows that the heavy quark propagátor acquires a mass AMq asymptotically. 
Thus 

< T 1V lV j '' ''''' 1 >~ _ e -AM Q (2T + 2L) (22) 

Large Wilson loops obey a perimeter law with o> = AMq. By restricting the discussion to 
classical fields the nasty issue of the divergences in the Coulomb energy does not arise. It 
would be interesting to find out, how ( p2]) compares with unquenched lattice calculations 
before and after cooling. A similar behavior is also present in the Schwinger model. 

In QCD, large Wilson loops are expected to obey a perimeter law. Does this mean 
that dilute instanton systems reflect faithfully on unquenched QCD ? We do not think 
so. In the presence of light quarks, the instantons and antiinstantons in the vacuum are 
screened over distance scales on the order of 1/2 fm. The screening is due to feedback of 
the quarks on the instanton configurations. While the effect is naively down by 1/N C , since 
the instanton density grows with N c , the screening persists to order AT?. This screening 
causes the topological susceptibility to be small and the rf to be massive. Dilute instanton 
systems in the presence of light quarks behave as a free gas. The operating mechanism 
is screening and not confinement. The lack of confmement is dramatically seen in (|8]) 
where the first term in the inverse correlator reflects on a dressed heavy-light quark cut. 
This problém occurs in all correlators. 

4. The above notwithstanding, the inverse correlator (|i3|) allows for an immediate 
interpretation in terms of effective interactions. In the long wavelength limit, the instanton 
size is small, and a local effective interaction between the effective fields Q and q (as 
opposed to the "fundamental" fields Q and q) can be derived much like the 't Hooft 
interaction between the light effective fields q. To show this, let (x, y) and (x', y') be 
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the coordinates of the light and heavy quarks respectively before and after they have 
interacted with an instanton (antiinstanton). From Bethe-Salpeter equation we read the 
vertex 

F(x, y, x', y') = -inN c J d A z JdU((x\[L] \x') ®{y\[H] \y') + 1-7) (23) 

This vertex function gives rise to an effective action Si. Averaging over the instanton 
orientation and position yields in the long wavelength (local) approximation 

The first bracket in (^4|) arises from the heavy quark part and the second bracket from 
the light quark part. Wick-rotating to Minkowski space gives 

which is to be compared with the 't Hooft vertex for two light flavors q = (u, d) 

/AM 2 \ 

£qq = {^J^J ( det + det ^qfl) ( 26 ) 

The ratio of the strengths is n q Q/n qq = AMq/2AM g ~ 0.08 and of order AT?, with 
«íQ ~ y/ň/N e . 

Similar arguments may be applied to QQ mesons as well. The induced effective La- 
grangian for heavy-light and heavy-heavy mesons is given by 

£ = q (i? - A M 3 ) q + (id t - AM Q ) Q 

/_l+7° l_l + 7° \ 

- (Q^-Qqq+^Q^-A^QqA^qJ 

í— 1 + 7° —1+7° l_l + 7° —1 + 7° \ 

- k qq Q-^-Q Q^-Q + iQ-^- Aa Q Q-^-A a Q (27) 
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where kqq = AMq/2nN c . For heavy baryons the construction applies as well. For qqQ 

systems the answer in Minkowski space reads 

(AM Q AM 2 \( —1+7° . 
£ qq Q = —y — 2n?řP~ I ' Q (detq L q fi + detq^q^ ) + 



\ Q^^A a Q (detq L A a qi? + detq fí A a q L ) ) (28) 



We note that n qq Q ~ n°/N 2 . A comparison with the conventional 't Hooft determinantal 
interaction shows that K qq Q/K qqq = K q Q/K qq ~ 0.08. The ratio of the strengths in heavy 
to light baryons is the same as the ratio of the strengths in heavy to light mesons. We 
note that both (p7|) and fl28|) are invariant under a spin fiip of the heavy quark SU(2)q 
and Ujy(l) violating, and that ( 28|) is chiral SU(2)l x SU(2)r invariant. 



In terms of the effective quark fields Q and q the dynamics is the one of the constituent 
quark model to order TOq. For heavy light mesons, the spectrum is 

M qQ = (m Q + AM Q ) + (m q + AM q ) + a qQ K qQ (l + ■ A%) (29) 

with a g Q = |0 9 q(O)| 2 ~ 1/OgQ, where 9 q(O) is the gQ-wavefunction at the origin and a q Q 
its typical size. Using the Van-Royen-Weisskopf construction JT2], we obtain 

N 

°«« ~ ^ff Q m 

with a q Q ~ 0.26 fm for fo < 290MeV ||. The spectrum (|2~9"D follows readily from 



the large distance asymptotics of the heavy-light correlator ( [12]) if we were to use the 
induced vertices (|25|) for simplicity. The interaction part in ( ^9|) corresponds to the scalar 
(1) and the (instanton-induced) Coulomb-like-exchange (A • A). The latter is leading 
and attractive in singlet configurations. Indeed, since the Casimir in the fundamental 
representation {Cp = A • A/4) scales with N c , the shift in the heavy-light mass is of order 
iV c ° with a qQ ~ iV c ° and K qQ ~ 1/N C , 

AM qQ C F a qQ K qQ ~ ^ OiqQKqQ (31) 



Recall that fl2"6|) in the large N c limit gives the following binding energy 

AM qq ~ -(1 - 4 S q ■ S q ) C F a qq K qq (32) 

For spin zero mesons 1 — 4S q • S q = 4 so that AM qq ~ —2N c a qq K qq . With this in mind, 
and using (|3~0D, the ratio of the bindings in heavy-light to light-light systems reduces to 

AM qQ l{AM Q ^(f qQ Ý (33) 



AM qq 8 V AM, ) \ fqq 
and similarly for the ratio of the bindings for heavy-heavy to light-light systems 

AMqq 1 / AMq \ 2 / f QQ ' 



mm 



AM qq 8\AM q J \f qq 

For D-mesons, we have AMd ~ 0.22 AM n (Jd ~ 290 MeV), while for charmonium 
AMj ~ 0.10 AM n (fj ~ 410 MeV). The drop in the interaction strength for heavy-heavy 
systems is compensated by an increase in the decay constant, in comparison to light-light 
systems. Since the light constituent mass is about 420 MeV, the pion binding energy 
is about 700 MeV. Thus the binding energy in heavy-light and heavy-heavy systems is 
about a 100 MeV. These estimates are consistent with the Coulomb estimates in the 



non-relativistic quark model [13 



Similar relations and estimates hold for baryons. They will be discussed elsewhere 
along with hyperfine splittings. Subleading effects in 1/N C arise from the left out parts in 
the color averaging, the non-planar contributions and the perturbative gluon exchanges. 
They are usually harder to assess without further assumptions. 

5. We have explicitly shown how the heavy meson correlator can be analyzed in powers 
of 1/rriQ prior to radiative corrections. Using a random and dilute instanton gas, we have 
shown that the heavy-light correlator can be calculated in the planar approximation and 
to leading order in the density fairly accurately. In the random gas the heavy and light 
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quarks acquire a mass to order mg. The induced mass of the heavy quark is related to 
the perimeter law displayed by the large Wilson loops. The latter reflects on screening 
in the gas. We have shown that in the long wavelength limit, the interaction induced by 
instantons between heavy and light quarks is attractive, and both chiral and heavy-quark 
symmetric. Our construction is immediately amenable to heavy baryons. It also allows 
for systematic investigations of the l/mg corrections and a possible assessment of the 
1/N C corrections. Some of these issues will be taken up in the future. 
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